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We analyze incommensurate charge-density-wave (CDW) and pair-density-wave (PDW) orders 
with transferred momenta (±Q,0)/(0, zbQ) in underdoped cuprates within the spin-fermion model. 

Both orders appear due to exchange of spin fluctuations before magnetic order develops. We argue 
that the ordered state with the lowest energy has non-zero CDW and PDW components with the 
same momentum. Such a state breaks C 4 lattice rotational symmetry, time-reversal symmetry, and 
mirror symmetries. We argue that the feedback from CDW/PDW order on fermionic dispersion is 
consistent with ARPES data. We discuss the interplay between the CDW/PDW order and d^ 2 _y 2 
superconductivity and make specific predictions for experiments. 


Introduction. The search for competitors to da^ 2_^2 
superconductivity (d-SC) in underdoped cuprates has 
gained strength over the last few years due to mount¬ 
ing experimental evidence that some form of electronic 
charge order spontaneously emerges below a certain dop¬ 
ing and competes with d-SC (Refs. [M) The two 
most frequently discussed candidates for electronic order 
are incommensurate charge density-wave (CDW) order 
(Refs, [ilill) and incommensurate pair-density-wave 
order (PDW), which is a SC order with a finite Cooper 


pair momentum Q (Refs. [32 


). Other potential can¬ 


didates are loop current orde r 139^ and CDW order with 
momentum near (tt, tt) (Ref. [40|). 

CDW order in underdoped cuprates has been proposed 
some time ago 0 and has been analyzed in detail by 
several groups in the last few years within the spin- 
fluctuation formalism [H 0, [ilhl IMil and within 
t — J model Emil. The initial discussion was focused 
on near-equivalence between d-SC and d-wave charge 
bond order (BO) with momenta (Q,Q) along zone di¬ 
agonal [ 1 ^, HO, but charge order of this type has 
not been observed in the experiments. It was later 
found [H, [ 23 , 0 , that the same magnetic model also 
displays a CDW order with momenta (Q,0) or (0,(5), 
which is consistent with the ran^ of CDW wave vec¬ 
tors extracted from experiments MSEliil. Such 
CDW order is also consistent with experiments that de¬ 
tect the breaking of discrete rotational and time-reversal 
symmetries in a (T, x) range where competing order de¬ 
velops In particular, when spin-fermion cou¬ 

pling is strong enough, the CDW order develops in the 
form of a stripe and breaks C 4 lattice rotational sym¬ 
metry. A stripe CDW order with ((5,0)/(0, (5) in turn 
gives rise to modulations in both charge density and 
charge current and breaks time-reversal and mirror sym¬ 
metries 23, 24, 28, 3 q| . 

The agreement with the data is encouraging, but two 
fundamental issues with CDW order remain. First, 


within the mean-field approximation, Tcdw is smaller 
than the superconducting Tc (and also the onset tem¬ 
perature for ((5, Q) order. It has been conjectured that 
Tcdw niay be enhanced by adding e.g., phonons 0, 
or nearest-neighbor Coulomb interaction or assum¬ 
ing the CDW emerges from already pre-existing pseudo¬ 
gap 0, 0. Tc,^is also enhanced by fluctuations be¬ 
yond mean-field [23|, [2J , but whether such enhancements 
are strong enough to make Tcdw larger than Tc remains 
to be seen. Second, stripe CDW order cannot explain 
qualitative features of the ARPES data away from zone 
boundaries 1361]. 


It has been argued [36| that ARPES experiments 
for all momentum cuts can be explained by assuming 
that the competing order is PDW rather than CDW. 
PDW order was initially analyzed for doped Mott insu¬ 
lators I 33 I. 371 . l38| , but it also emerges in the spin-fermion 
model [28[ with the same momentum ((5,0)/(0, (5) as 
CDW order and its onset temperature Tpdw is close to 
Tcdw (the two become equivalent if one neglects the cur¬ 
vature of fermionic dispersion at hot spots 27, 0). 
Given that PDW order explains ARPES experiments, it 
seems logical to consider it as a candidate for competing 
order. Just like CDW, the PDW order develops in the 
form of a stripe and breaks C 4 lattice rotational sym¬ 
metry 28, 3^, if, again, the coupling is strong enough. 
However, it does not naturally break time-reversal and 
mirror symmetries [35| (although it does so for a par¬ 
ticular Fermi surface geometry [3^), and the mean-field 
Tpdw is also smaller than Tc for d-SC. 


In this communication we build on the results of the 
generic Ginzburg-Landau analysis 0 and propose how 
to resolve the partial disagreement with experiments for 
pure CDW or PDW orders. We first re-iterate that pure 
CDW/PDW orders emerge in the forms of stripes only if 
the spin-fermion interaction g is strong enough. In prac¬ 
tice, g has to be at least comparable to the upper energy 
cutoff of the spin-fermion model A (see details below). 
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FIG. 1. The Brillouin zone, the Fermi surface, and the hot 
spots. We label bonds connecting hot spots as A, B, C, D, 
a and b. Inset: the structure of the mixed CDW/PDW state 
in one of the hot regions. 


For smaller couplings the system develops a checkerboard 
order for which C 4 symmetry is preserved [i^ . The spin- 
fermion model is a low-energy model and it is rigorously 
defined only when the coupling g is smaller than A. In 
this respect, stripe CDW or PDW orders emerge, only 
at the edge of the applicability of the model. Here we 
consider spin-fermion model at smaller couplings, well 
within its applicability range, and allow both CDW and 
PDW orders to develop. We show that the system devel¬ 
ops a mixed CDW/PDW order, in which a CDW com¬ 
ponent develops between hot fermions separated along, 
say, Y direction and a PDW component develops between 
fermions separated along X direction (see Fig. 1). Be¬ 
cause the momentum carried by an order parameter is 
the transferred momentum for CDW and the total mo¬ 
mentum for PDW, the CDW order along Y and the PDW 
order along X actually carry the same momentum (0, Q). 
We argue that such a state further lowers its Free en¬ 
ergy by developing (via an emerging triple cou plin g) sec¬ 
ondary homogeneous superconducting orders [28|. This 
effect favors the mixed CDW/PDW state over the pure 
checkerboard CDW or PDW states, which would oth¬ 
erwise all be degenerate. The mixed CDW/PDW state 
breaks C 4 symmetry because both orders carry either 
momentum (Q,0) or (0,(5), but not both, and it also 
breaks time-reversal and mirror symmetries as the pure 
stripe CDW order with ((5,0) or (0,(5) does. 

The presence of PDW component is relevant for the 
interpretation of the ARPES data. Without it, the 
fermionic spectrum in the CDW phase would contain 
the lower energy branch, which never crosses Fermi level, 
and the upper energy branch, which would approach the 


Fermi level from above as the momentum cuts enter the 
arc region. As discussed in [s^, this is inconsistent with 
the data @ which show that the dispersion approaches 
the Fermi level from below. We show that the presence of 
PDW component changes the structure of fermionic dis¬ 
persion in such a way that now the lower branch crosses 
the Fermi level in the arc region (see Fig. 2), in full agree¬ 
ment with ARPES experiments. 

We also consider the interplay between CDW/PDW 
order and d-SC and present the phase diagram in Eig. 3. 
The reduction of the superconducting Tc in the coexis¬ 
tence region with CDW/PDW is the obvious consequence 
of competition for the Eermi surface. A small (of order 
g/A.) drop of Tc upon entering the coexistence region is 
the result of a weak first-order CDW/PDW transition. 
There exists, however, a more subtle feature of the phase 
diagram. Namely, a secondary SC order is generated by 
CDW/PDW order, which preserves the same sign of the 
gap along each quadrant of the Eermi surface. Below 
Tc for d-SC, this secondary superconducting order cou¬ 
ples with dx 2 _y 2 order, and the net result is the removal 
or shifting of the gap nodes. Simultaneously, the CDW 
order acquires an extra component with s-form factor, 
i.e., the magnitude of its s-wave portion increases. We 
propose to verify these through experiments. 

The model We follow previous works [SHIillii 
and consider emerging charge order within the spin- 
fermion model [ 4 ^. This model describes interactions 
between itinerant electrons and their near-critical an¬ 
tiferromagnetic collective spin excitations in two spa¬ 
tial dimensions. Eight “hot” spots, defined as points 
on the Eermi surface separated by antiferromagnetic or¬ 
dering momentum (7r,7r) (points 1-8 in Eig. 1), are the 
most relevant for destruction of a normal Eermi liq¬ 
uid state. The known instabilities of the s pin -fermion 
model include d-SC (e.g. (cice), see Eig. 1) fiol. Id^. . 
bond charge order (BO) with momenta (±(5, i(5) (^-g- 


(cbe)) [ii,0l23, CDW order with momenta (0, ±(5) 
and (±(5,0) (e.g. {c\c 2 )) [23,[2l,[3l| and PDW order with 
momenta (0, ±(5) and (±(5,0) (e.g. (C 1 C 2 )) [ 23 , 0 . The 
model has an approximate SU(2) particle-hole symme- 
try [ 19 I, [ 2 ^, [ 27 I, 0, [ 13 , which becomes exact once one 
linearizes the fermionic dispersion in the vicinity of the 
hot spots. This gives rise to near-degeneracy between 
d-SC and BO and between CDW and PDW. 

The Ginzburg-Landau analysis We introduce four 
order parameters: for SC, 4> for BO, for PDW, and 
p for CDW respectively. SC and BO order parameters 
connects hot spots along diagonal bonds, which we label 
as a and 6, while PDW and CDW connect hot spots along 
vertical and horizontal bonds, which we label as A-D in 
Eig. 1. We define the CDW order parameter on bond 
A as Pa {c\^ 2 ) and use analogous notations for other 
order parameters. The effective action has three terms: 


‘^cdw/pdw[P? ± ‘^sc/bo[^: ± ‘^int (1) 
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FIG. 2. Fermionic dispersion in the antinodal region in the 
presence of the mixed CDW/PDW order. Upper panel - the 
dispersion in the presence of CDW/PDW order for various 
kx {kx = TT corresponds to the cut along the Brillouin zone 
boundary). Middle panel- the spectral function. Thin line 
on both panels is the bare dispersion. Bottom panel - exper¬ 
imental data from Ref. for comparison. The experimental 
data have been taken below Tc and show a gapped dispersion 
in a wider range of tv — kx. 


The Scdw/pdw[p^ term is of our primary interest. Keep¬ 
ing the SU(2) symmetry exact, we follow Ref. [1^ and 
combine PDW and CDW orders on a given bond (say, 
bond A) into a 2 x 2 matrix order parameter 

( 2 ) 

where pA ^ c|c 2 , 'i^a ^ ciC 2 , and Ua is a SU(2) 
matrix “phase”. The order parameters Ab,c,d and 
phases Ub,c,d are similarly defined (see Supple¬ 
mentary Material (SM) for details). Minimizing the 
Free energy, we obtain F = Tt{UaUIjUbUd) = 
-2 , + = a/IpbP + IV^bP = |A^|, and 

y/|pcp + = |Aa:|. Under these 

conditions, the CDW/PDW action becomes 

5cdw/pdw =f (|A.|^ + |A,|^) + /3(|A.r + |A,r) 

+ (^-^)|A.p|A,p + 0(A^) (3) 


where Of ^ Ajvp x {T ^cdw)/^cdw and T/dw — ^pdw ^ 9 
(Ref. [i^). The prefactors and p are determined by 
different convolutions of four fermionic propagators (the 
square diagrams [H, HH). At g A we have p ^ 
l/{i;|A), j3 - \og{K/g)/{v%K), and j3 ~ {A/g)/{v%A). 
We see that p is the largest term, hence the action (j3j) is 
minimized when |A| = |Aa,| = |A^|. Because P — p <f)^ 
the action is unbounded, which implies that the tran¬ 
sition is first-order and sixth-order terms (coming from 
six-leg diagrams) have to be included to stabilize the or¬ 
der. Including these terms we obtain a first order into 
CDW/PDW state at Tcdw/pdw = Tcdw(l + 0{g/A)). We 
emphasize that this temperature is higher than the one 
for a pure CDW (or PDW) transition. 


The constraint F = aUIjUbUd) = —2 leaves the 

ground state hugely degenerate - the order parameter 
manifold is SO(4) x SO(4) (Ref. [H). This manifold in¬ 
cludes pure CDW and pure PDW checkerboard states 
and mixed CDW/PDW states. To select the actual 
ground state configuration we note that, if CDW and 
PDW orders have components which carry the same mo¬ 
mentum Q, the Free energy is further lowered by creating 
a secondary order whose magnitude is a product of CDW 
and PDW order parameters. This secondary order is a 
homogeneous SC with equal sign of the gap along each 
quadrant of the FS [1^ One can straightforwardly check 
that the reduction of the Free energy is maximal when in 
a nominally checkerboard state CDW occurs along ver¬ 
tical bonds and PDW occurs along horizontal bonds or 
vise versa, i.e., each order develops in the form of a stripe. 
This corresponds to either P’a^b = Pc,d = 0 (as in the 
inset of Fig. 1) or p)c,D = Pa,b = 0, the choice breaks C 4 
lattice rotation symmetry. Furthermore, the stripe CDW 
order parameters pA and pB and PDW order parameters 
p)c and 'ijjD get separately coupled by fermions away from 
hot spots, and the coupling between pA and pB locks the 
relative phase of pA and pB such that pB = AipA (Ref. 
M)- The choice of the sign breaks time-reversal and 
mirror symmetries. The coupling between 2 pc and p)D 
does not lock their phases. 

Feedback from CDW/PDW order on fermions We 
now show that the feedback from stripe CDW/stripe 
PDW order on the fermionic dispersion At k ^ ('7r,0), 
taken as a function of ky for various kx = tt — 6 kx, 
yields results in quite reasonable agreement with ARPES 
data [1, [ 1 ^ . Previous studies have shown [ 2 ^ that a pure 
CDW order can explain the ARPES spectrum for a cut 
along the BZ boundary, but not for cuts that are closer 
toward BZ center (see Ref. 3^,113). To obtain the disper¬ 
sion along various cuts in the presence of both CDW and 
PDW, we have extended our analysis of the CDW/PDW 
order to a finite momentum range away from the hot 
spots. We find that at the BZ boundary, the CDW order 
has a larger amplitude due to better ES nesting but the 
PDW component increases as the cuts move towards the 
hot spots. We present the details in SM and show the 
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FIG. 3. The phase diagram. The transition into CDW/PDW 
state is weakly first-order and the superconducting Tc drops 
by a finite amount upon entering into coexistence region. In 
the region labeled as “pre-emptive” discrete Ca and time- 
reversal/mirror symmetries are broken but continuous U{1) 
translational symmetry (associated with the locking of the 
common phases of pA and pb) remains unbroken [^. In the 
shaded region, Mott physics develops and the onset temper¬ 
ature of charge ordering shrinks. 


results in Fig. 2. There are three key features in our sce¬ 
nario that are qualitatively consistent with experiment: 
(1) at the BZ boundary {kx = tt), the locus of mini¬ 
mum excitation energy shifts from kp to a larger value 
kc ~ Q/2, where Q is the CDW momentum, (2) as kx de¬ 
creases, the excitation approaches the Fermi level from 
below, and (3) at kx when the Fermi arc emerges, the 
fermionic dispersion becomes flat for \ky\ > kp- These 
features are also reproduced by pure PDW order and 
from a spatially homogeneous self-e nerg y arising from a 
d-wave CDW order peaked at (tt, tt) [21|. However, both 
these scenarios do not immediately explain the obser¬ 
vation of broken time-reversal symmetry or CDW order 
with small incommensurate momentum. To obtain quan¬ 
titative agreement with the experiments, we would need 
to know how CDW and PDW order parameters depend 
on frequency. This would require one to model the bare 
dispersion far away from kp and solve complex integral 
equations for frequency-dependent order parameters. 


Interplay between CDW/PDW order anddx‘2-y2 super- 
eonduetivity We next consider other terms in the effec¬ 
tive action in Eq. ([1]). The term S^c/ho das been analyzed 
in [23,[23,[3 i|. When SU(2) symmetry is exact, d-SC and 
BO orders are degenerate and the action has four Gold- 
stone modes. Once SU(2) symmetry is broken by FS cur¬ 
vature, only d-SC order develops below Tc. We assume 
that this is the case and keep only d-SC component 4^ 
in 5sc/bo, i-e- reduce it to 5sc/bo = with 

Os ~ X (T - Tc)/Tc, Tc ~ g, and ps ~ A/(uFff)^- 


The coupling between CDW/PDW and d-SC orders is 
again obtained by evaluating the square diagrams. The 
calculation yields Si^t = with ^ l/{v‘/g). 

Note that the magnitude of the coupling is phase sensi¬ 
tive, hence the phase locking between pA and at ±7r/2 
is important (see SM for details). 

The analysis of the full action is straightforward and 
we show the results in Fig. 3. The mean-field tempera¬ 
ture Tcdw/pdw ^ ^cdw is comparable to Tc near the SDW 
boundary but is enhanced by fluctuations 23,0, [ 2 ^. We 
assume that this enhancement lifts Tcdw/pdw above Tc at 
large Because CDW/PDW transition is first-order, Tc 
jumps upon entering into the coexistence region, but the 
jump is again small in g/K. Similar behavior has been 
recently observed in Fe-pnictides At small T, the 

CDW/PDW and d-SC orders coexist. 

The phase diagram in Fig. 3 is similar to that for pure 
CDW order [i^], but there are some extra features. First, 
the combination of CDW/PDW orders induces a sec¬ 
ondary SC order [ 2 ^ with a non-zero gap along zone 
diagonal ( 5 -wave or dxy). In the coexistence region the 
order 4^^ couples with d-SC order 4^ and, as a result, 
gap nodes either get shifted {d s state) or removed 
{d + e^^5 state). A similar coupling has been examined 
in the context of the Fe-pnictides [52[ . A finite gap along 
zone diagonals has been observed in ARPES at doping 
X < 0.1 (Ref. [13) and also inferred from Raman spec¬ 
troscopy [HI]. Second, by the same logic, the d-SC and 
PDW orders induce a secondary 5-wave CDW order with 
the same momentum as the primary one. We propose to 
search for SC gap opening or node shifting and to ex¬ 
amine the 5-component of CDW order in the coexistence 
region. 

Conelusions In this letter we proposed a state with 
unidirectional CDW and PDW orders which carry the 
same momentum. We argued that this state is a mem¬ 
ber of the ground state manifold of the low-energy spin- 
fermion model and its energy is further reduced by in¬ 
duction of a secondary SC order. We further argued 
that CDW/PDW state has a number of features con¬ 
sistent with experiments: it breaks both C 4 and time- 
reversal symmetry and the feedback from CDW/PDW 
order on fermions reproduces the ARPES data from the 
BZ boundary to the tip of the Eermi arc. The transition 
into CDW/PDW state is weakly first-order and occurs at 
a higher transition temperature than that for a pure uni¬ 
directional CDW or PDW orders. We considered the in¬ 
terplay between CDW/PDW order and d-SC, and found 
that a SC gap becomes non-zero along zone diagonals. 
We proposed to search for this gap opening in the region 
where charge order and d-SC coexist. 

We acknowledge useful discussions with W. A. Atkin¬ 
son, D. Chowdhury, E. Eradkin, S. Kivelson, P. A. Lee, 
C. Pepin, and S. Sachdev, The work was supported by 
the DOE grant DE-EG02-ER46900 (AC and YW) and 
by NSE grant No. DMR-1335215 (DEA). 
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Supplementary Material 


I. DETAILS OF THE GINZBURG-LANDAU ACTION 


A. CDW/PDW sector 


The Ginzburg-Landau (GL) action for the CDW/PDW order has been derived and studied in detail in Ref. [l|. We 
briefly review the analysis here and apply it for our purposes. 

When the curvature of the Fermi surface (FS) at hot spots can be neglected, spin-fermion model has SU(2) symmetry 
which makes CDW and PDW orders degenerate. The action of the spin-fermion model can be rewritten in an explicitly 
SU(2)-symmetric form in terms of particle-hole doublets at each of the hot spots 1-8, 


=(l(-t)' 

=( %k)) ■ 


C2t(fc) ) 
c\^{-k) \ 
C4t(fc) ) 

C6t(fc) \ 

44.(~^) / 

C8t(fc) ^ 

44“^) / ’ 


( 4 ) 


In this notation, the CDW order parameter p’s and the PDW order parameter t/’s involving the same pair of hot spots 
(labeled as A, 5, C, and D in Fig. 1 in the main text) can be combined into a 2x2 matrix that couples bilinearly to 
particle-hole doublet The four 2x2 matrices are 


Ar = 


Pa 

-PA IpA 




V’B PB 
~P*B '^B 




'^C Pc A ( ‘4’D Pd\ 

-PC rc)' ^ V -Pd ro ) ■ 


( 5 ) 


For convenience, we also define SU(2) phases Ua,b,g,d via Ua,b,c,d = ^a,b,c,d/\/\pa,b,c,d\^ + |V’a,b,c,bP- Each 
order parameter changes sign under a momentum shift of (tt, tt) (e.g., between the pair 1,2 and the pair 3,4 in Fig. 
1 in the main text) because spin-mediated interaction is repulsive. The magnitudes of the CDW and PDW order 
parameters between 1,2 and 3,4 do not have to match as these two pairs of hot spots are not equivalent. For simplicity, 
below we neglect this non-equivalence and assume that order parameters just change sign under a momentum shift 
by (tt, tt) (this is often termed the d— wave approximation for the form-factor of the charge order). 

The full effective action in terms of CDW and PDW order parameters up to quartic order is 


5 eff =a' Ti{A\Aa + A^Ab + A^Ac + A^j^Ad) 


- (/i + h) Tr(AAA]4A^ A^ + AbA^AbA^ + AcA^AcAI + AbA^AbAI,) 


- 2/3 Tr 


'(A^A], + AbA^)(AcA^ + AbA|,)] - 4/4 Tr (aJ^AbA^Ab) 


=2a'i\pAf + + \pBf + li^Bf + \pcf + IV’cl' + \PDf + I^dI") 

- 2(Ji + I 2 ) [(IpaP + + (I/Bp + IV’Bp)^ + (l/Cp + IV’Cp)^ + (I/Bp + If^Bp)^] 

- 4J3(|p^P + I^Ap + I/Bp + |^Bp)(|/Cp + IV’Cp + l/Bp + IV^Bp) 


- 4/4\/(|paP + |fi’AP)(|/BP + |fi’BP)(|/cP + |fi’cP)(|/BP + \^D\^)Tk (u^UcUiUD) 


(6) 


In mean-field analysis the coefficient a' is proportional to g~^ — n(T), where ^eff is the effective four-fermion 
interaction and 11 is the polarization operator (see Ref. 0 for details). The polarization bubble 11 (T) increases as 
temperature decreases, and a' = a'{T) changes sign at the CDW instability temperature T = Tcdw = Fpdw- By 
dimensional argument 11 ^ A/vpf{T/Tcdw)^ where A is the upper cutoff of the spin-fermion model, and f{x) is a 
dimensionless function. Then a' = A/4(/(l) -/(T/T^dw)) = A/4[(T - Tcdw)/Tcdw]/'(1) - A/4[(T - Tcdw)/Tcdw]. 

The calculation of Tcdw = Tpdw requires more care. If one neglects momentum and frequency dependence of 
one obtains that Tcdw is non-zero only if ^feff exceeds some critical value iii- This is a consequence of the velocities 
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FIG. 4. The diagrams for coefficients Ii - I4. 


of hot fermions at and ki? + Q being generally not antiparallel. However, once one includes the fact that spin- 
fermion interaction ^eff = 9 X{q^^) is mediated by a boson with near-divergent dynamical susceptibility x(^,co’), one 
obtains Q that the threshold vanishes when the magnetic correlation length diverges. In this limit, the CDW/PDW 
instability occurs for arbitrary values of the spin-fermion coupling and Tcdw = Ppdw ^ 9 - The spin-fermion model is 
justified as a low-energy model when interactions do not take a fermion outside of the low-energy subset, which holds 
when g A. 

The coefficients /i, 2 , 3,4 are obtained by evaluating the four square diagrams in Fig. 01 In explicit form, 

LOm 

/ 2 =-T^ / ^Gl{uj^,k)Gl{uJm,k), 

UJm 

h = -Tj2j ■^Gl{LUm,k)G2{com,k)G5{com,k), 

UJm 

id = — T J -^^Gi{(Jrn^ k)G 2 {uJrn^ k)G 5 {iUrn^ k)GQ{iUrn^ (7) 


where Gi{uJm^ k) is the Green functions for a fermion near hot spot i, and momentum k is defined as a deviation from 
this hot spot. The evaluation of the integrals yields im 


ii=- 

/2 = 0 , 


1 


1 


Wn'^vlvy Ak ’ 


r 1 1 , ^ 

I3 = o o — log — 


I 4 


IGn^vlvy Ak 

1 


32VxVy T 


(8) 


where Vx and Vy are x and y components of the fermi velocity Vi? at hot spot 1 , and is the momentum cutoff 
^ A/vp. Using the fact that T ^ we have Ii ^ l/(u|.A), Is ^ l/(u|.A) log(A/^) and I 4 ^ 

Because U\UcUq^d ^ SU(2), the term proportional to I 4 in Eq. ([ 6 ]) is minimized when Tt(^\UcUq^d^ = 

— 2 . We also find that in the ground state a/\pa\^ + = \/\pb\^ + and aJ\pc\^ + = 


a/\Pd\^ + = |Aa,|. Under these conditions, Eq. ([ 6 ]) becomes 


^cdw/pdw 


^ n X 

= -( A.^ 


\Ayf) + l3i\A.\^- 


\Ay\^) + {^-m,f\Ay\ 


0(A®), 


(9) 


where a = 4(a', [3 = 4|/i -h / 2 I, = I 6 I/ 3 I, and (3 = 8 |/ 4 |. This is Eq. (3) of the main text. 

As we said in the main text, P > P and the action is unconstrained at the quartic level. We assume that sixth 
order terms, given by convolutions of six Green’s functions have positive a prefactor 7 and constrain the action. An 
order of magnitude evaluation of 7 gave 7 ^ l/{vpT^) ^ l/{vpg^). 
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B. Stabilization of the mixed CDW/PDW order 

In an SU(2)-syninietric model, all states that satisfy |Aa,| = |Ay| = |A| and F = Tt{UaU^UbUI)) = —2 are 
degenerate ground states. These include pure checkerboard CDW and PDW states and a mixed CDW/PDW state 
with stripe CDW and PDW orders between pairs of hot spots along x and y directions respectively (or vise versa). 
We show that, when FS curvature is non-zero, the mixed CDW/PDW order generates a secondary homogeneous 
SC order, and this favors the mixed CDW/PDW order over pure CDW or PDW states. We also show that the FS 
curvature lifts the degeneracy between CDW and PDW orders and makes both quadratic and quartic terms in the 
effective action anisotropic. We show that this also favors the mixed “stripe” CDW/PDW state for some range of 
parameters. 


1 . Coupling to secondary homogeneous SC order 


As we said in the main text, in the mixed state, CDW and PDW orders which carry the same momentum generically 
induce a secondary homogeneous SC order via triple coupling terms, and these terms can lower the Free energy. We 
show that the mixed CDW/PDW optimizes such coupling. The induced secondary SC order was shown in Ref. [l| to 
be a mixture of s-wave and d^^^-wave. We define s-wave SC and d^^^-wave SC order parameters as and The 

action for the secondary SC order is given by [I| 


V’] =y + ^dJiPArn + Pd^a + PlVc + pars) 

+ - ^dAiPsi’D + poi’i + P*Ai’c + Pc^a)] + h.c. 

= + ^d^C + \^s - ^dj^) - + ^d^,){pA^*D + Poi’A + P*B^C + pci’s) 

+ A, - ^dJiPBrn + Pd^b + pWc + PcV’a)] + h.c., (10) 

where ( 0 ^ 5 )“^ > 0 is the susceptibility of the secondary SC orders (for simplicity, we take this susceptibility to be the 
same for and ^ qualitatively similar mixed CDW/PDW ground will still result if the two susceptibilities 

are not the same). When both p and ^ are nonzero, superconducting orders and are induced and the Free 
energy is lowered. 

To minimize the Free energy m, we maximize the magnitude of the two combinations of CDW/PDW order 
parameters pA’i’h + Pd^’a + PbV’c + Pci’s Psi’s + Poi’s + P*Ai’c + Pci’ a- We define 


PA = |A| cosPac^’^'’, 
PB = |A| COS^bC*'^®, 
pc = |A| cos6'ce*'^^, 
Pd = |A| cos 9 


i^A = |A|sin0^e'^^ 
i)B = |A| sin 6»Be®^®, 
i)c = |A| sin^ce*^*^, 
V’D = |A|sin 0 De*^®, 


( 11 ) 


where 0 < 9a,b,c,d < -tt < 4>a,b,c,d < tt, and —tt < (j)A,B,c,D < tt. 
The first CDW/PDW combination term in Eq. (pd)|) hence becomes 


PArD^PDrA+PBrc^pcrB 

= |Ap cos^A + cosdy) -t-cosd^ + cosdc sindBe*^"^^”^^^ 

= |Ap [cosd^ sind^e*"^^ + cosdy) sind^e*"^^ + cosd^ sindce*"^^ + cesdc sind^e*"^^] 


= |Af 


je^Wi 


+e 


- 02)/2 

— h03—04)/2 


sin(dy^ + Ob) cos 
sin(dc' + d^) cos 


01 + 02 
2 

03 + 04 


— i sin(dA — Ob) sin 
— i sin(dc — d^) sin 


01 + 02 
2 

03 + 04 


]} 


( 12 ) 


where in the third line we have defined 0i = 0 a “ 0d , 02 = 0 d + 0a , 03 = 0 b + 0 c ? and 04 = 0c “ 
For generic 0i,2,3,4 the magnitude of Eq. m is maximized when 

Oa = 0, Ob ~ "2 ’ A ^ D 

9b = 0 , 0 c = ^, or B^C. 


(13) 
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Repeating the same arguments for the second combination term + Pd'^b + Pa^c + Pc'^A obtain one more 

set of conditions 

Ob =0, Od = ^, or AD 

0a = 0, 0c = ^, or B^C. (14) 

Combining, we obtain 

8a=0b= 0, ec = eD = ^, or A, B^C,D. (15) 


Then pAiph + + Pg^b = 1^1' 


A{(f>A-4>D) 


+ e 




and psro + Pd^b + pWc + P*c^l = 




. Maximizing the magnitude of both these terms, we obtain the conditions on (p^s and 


0’s as 


(pA + (pB + (pc - (pD = 0. 


(16) 


It is easy to verify that Eqs. and yields \pa\ = \pb\ = l^cl = 10^1 = |A|, pApB'ipcP^h = or 

I Pc I = I Pd I = 10^1 = 11 = |A|, pApBp^ci^D ~ \paPbP^cP^d\- These states are exactly the mixed CDW/PDW state 
we described in the main text, related by a lattice C 4 rotation. 

Before we proceed, we remind that P’s, 0’s, and 0’s are not free parameters - they are constrained by the condition 
r = TtCUaU^UbUI^) = —2, and Eqs. and have to be consistent with this condition. We recall that in our 
notations 


Ua 


Uc 


sinP^e^^^ cosP^e-^^^ 

— cos Oag'^^^ smOAe~'^^^ 

sinPce*^*^ cos 

— cos sin 0 ce~^^^ 


, Ub = 


, Ud = 


smUBe^ cost 
— cos0bc~'^^^ s\a0bc~^^^ 


sin cost 

- cos sm. 0 DC~'^^^ 


(17) 


Plugging Eqs. and into Eq. ([T7j) we find that indeed P = Tt{UaU^UbUI)) = —2. Hence, the mixed 
CDW/PDW state is truly a ground state. 

So far for the minimization with respect to P’s we have assumed that 0i,2,3,4 are completely generic. Eor the special 
case when 0i = —02 and 03 = —04, the condition on P’s are less strict - from Eq. m we see that in this case one 
only needs to satisfy 


Oa = 7r/2 — P^, and Ob = 7r/2 — Oc = P- 
01 = —02 = 03 = —(pA‘ 

Doing the same for PbP^d + PdP^b + Pa^^c T Pc'P^a and combining with ([18]) we find that 

Oa = Ob = 7r/2 — Oc = 7r/2 — Ob = 0^ 

(pA — = -(pD - 0A = -(pB - = (pc - 0D, 

(pB — = (pD — ^B = —(pA — 0c = —(pc — 0A- 


(18) 

(19) 

( 20 ) 


The configuration in Eqs. m and dsni) would give the same ground state energy as our mixed CDW/PDW state 
and hence has to be considered. However, again one needs to check if Eq. ([20|) is consistent with P = — 2, which in 
terms of P, 0’s, and 0’s becomes 


Re 


- sin^ 20 
4 


gd0C-0A)^ ^p(f>B-(f>D)^ 


— ^sin^ Oe 40a+0c) _ gQg 2 ti(pA-\-(pc) 
Eor a generic P, this condition is satisfied if 




2 Q^iiPB-C^D) _ gog2 Q^i{4>B-ci^B)^ 


- 1 . 


( 21 ) 


<pA + (pc - <pA - (pc =7r, 

^B — (pD — + (pB =7r, 

^A-\- (pc — ^B (pD =0. 


( 22 ) 
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or 


(t>A + (t>c - (t>A- 4>c =0, 

— 4>d — 4^b = 0 - 

^A + (I^C - ^B + (t>D =0. 

- ^A = 7 r, 

— ^B = 0 . 


(23) 


We have verified that neither Eq. nor ([231) consistent with Eq. (EU]) . Hence Eqs. ([TU]) and (^0]) do not 
correspond to a true ground state of the CDW/PDW order, and only Eqs. ([T^ and (IT^ are the correct conditions 
that minimize the Eree energy in Eq. m- Minimization of Eq. ([TQ|) with respect to and yields |^s| = 

EjAp/a^, and <Stri = —SI3\A\^^ where ^ Y‘^ ja^. We can then effectively write down the Free energy m 
as 


Stri = —^^|A|^(cos6>a cos Ob sinOc smOo + sihOa smOs cos Oc cos Od)- 


(24) 


2. Other effects beyond hot spot approximation 


Two other effects beyond hot spot approximation have been considered in Ref. [l|. First, at the quadratic level, 
there appears the term ^ — 6 a\ 2 p\‘^ in the action which favors PDW order. Second, at the quartic level, coupling to 
fermions away from hot spots reduces the CDW Free energy by ^ —^ 7 \pa\‘^\pb\‘^ and favors CDW order. These two 
effects and the one considered in the previous subsection give rise to an extra piece in the effective action of the form 

6S = — fe|A|^(sin^ Oa + sin^ Ob + sin^ Oc + sin^ Od) 

— 5P\A\^{cos0a cos Ob sinOc smOo + sihOa sinOB cosOc cos Od) 

— ^7| A|^(cos^ Oa cos^ Ob + cos^ Oc cos^ Od), (25) 


where ^/3, ^7 > 0. In the ground state we have Oc = Od = 0^ and Oa = Ob = Oy. The values of Ox and Oy depend 
on the parameters 5a^ 5/3, and J 7 . If 5a term is the largest, the ground state has Ox = Oy = 7 rl 2 and is a pure PDW 
state; if 5j term is the largest one. Ox = Oy = 0 and the ground state has pure CDW. If 5f3 is the largest, 0x,y = 0 , 
Oy^x = 7 r /2 and the ground state is our mixed CDW/PDW. For this case it is easy to verify that the condition on 5/3 
to stabilize a mixed CDW/PDW state is 


5/3 > 


5^ 


5a 


(26) 


We note that when this condition is satisfied, the ground state will no longer have equal magnitudes of the CDW and 
PDW orders (however, the ground state will remain a mixed CDW/PDW state). It is relevant for this reasoning that 
the transition into the mixed CDW/PDW state is first-order because if it was second-order, the condition ([26|) could 
not be satisfied at temperatures right below Tcdw and the system would first develop a pure PDW order. 


C. SC/BO sector 

In this subsection we analyze the interplay between (ia, 2 _^ 2 -wave SC and bond charge orders (BO) with momenta 
{Q,3zQ). We define SC order parameters ^ (ciCe), T 5 ^ (C 2 C 5 ), and BO parameters 4>a ^ {c\cq), ^ (^ 2 ^ 5 ). All 

order parameters change sign when fermionic momenta are changed by (tt, tt), namely, (cice) — (csCg), (C 2 C 5 ) — (04 C 7 ), 
(cjce) — (C 3 C 8 ), and (C 2 C 5 ) — {c\cj). For SU( 2 )- symmetric model the effective action can be written as iii 

'5sc/bo =<^s (l^al^ + I^a7 + (l^'bP + l^bP) + /?« (l^'aP + I^a7 + /?« (I^'bl^ + l^’bP) • (27) 

This action has 0(4) x 0(4) symmetry. The two SC order parameters and can be combined into dx 2 -y 2 and 
B 2 g order parameters as 


= («'a + ^'b )/\/2 and = (^'a - ^b)/V 2 . 


(28) 
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Equivalently we have ^6 = {^d - 

By obvious practical reasons, we only consider dx‘ 2 -y 2 wave SC and set zero. We then have 

^2 = The action takes the form 


‘5sc/bo —<^s(|Tc^P + l^aP + l^6p) +/^s 

^ a .(| v|/,|2 + | i.„|2 + |$,| 2 ) + ^ (|,|>,|2 + |$^|2 ^ ^ ^ _ |^^| 2)2 _ 

The fist two terms in the last line of Eq. (E9|) describe a model with 0(6) symmetry The last term breaks the 
0(6) symmetry and gaps out would be longitudinal Goldstone mode between \^a\ and As a result, the model 
of Eq. ([29|) has four Goldstone modes instead of five for an 0(6)-symmetric model. Because /3s > 0, for the ground 
state of Eq. ([29|) we have \^a\ = |^6|- Defining ^ = \/2^a and setting \^a\ = \^b\^ we obtain from ([29|) 



5sc/bo[«'d, $] =aa(|4'd|' + + Y ■ 


(30) 


In terms of and the ground state has an 0(4) degeneracy. 

The degeneracy is lifted once one includes into consideration the ES curvature at hot spots. Then the action 
becomes 


>Ssc/bo[^d,$] =aa|^dp + 5,|$|2 + ^(|^,|4 + |$|2)", (31) 

where < d^. In this case the pure SC state minimizes the Eree energy and becomes the true ground state of 
the SC/BO sector HEi. BO always comes second and, because this order does not additionally break any discrete 
symmetry, there is no possibility to lift the instability temperature for this order above the superconducting Tc. 


D. The coupling between d^2_y2-w8L've SC and the mixed CDW/PDW order 

In this subsection, we discuss the interplay between da, 2 _^ 2 -wave SC and the mixed CDW/PDW order. We first 
show that at low temperatures, but in the range where GL expansion is applicable, the CDW/PDW gap A is 
parametrically larger than the SC gap T. To see this, we note that the GL action for mixed CDW/PDW is <Scdw/pdw = 
(a|Ap —/3|A|^+7|A|^. Minimizing this action we obtain |Ap = (^ + y//3^ — 3ja)/3^. Below Tcdw, is negative and 
scales as —A/v^. Then A^ ^ g^^JA/g. Meanwhile, the SC gap is given by |Tp = |(as|/(2/3s), where 

\as\ ^ A/vp and (3s ^ A/{vfTc)^ ^ A/{vFg)^- Using the fact that Tc and Tcdw are both of order g, where, we recall, 
g is spin-fermion coupling, we obtain |Tp ^ \A‘^\y^g/A |Ap. This result implies that the interaction with SC 
order does not distort substantially the inner structure of the CDW/PDW order. In particular, the condition T = —2 
yields the minimum of the Eree energy also in the presence of < 3 ^, 2_^2 SC. 

We now derive the effective action. We first note that da, 2_^2 SC and CDW/PDW orders do not couple via triple 
couplings ^ Tp?/, because all three order parameters change sign when fermionic momenta are shifted by (7r,7r), and 
the contributions from fermionic momenta that differ by (7r,7r) cancel out. Rather, da, 2_^2 SC couples separately to 
CDW and PDW components. We will drop the subscript d and use T for dx‘^-y‘^ SC order parameter. 

Eor PDW components, we find two types of couplings to SC, which are given by the diagrams shown in Eig. [SJa). 
The first diagram yields a coupling term — Is where 

C d‘^k 

I5=-TJ2 / -^G^iuJm,k)G2iuJrr^,k)Gl{-Um,-k). 


(32) 
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(b): CDW and SC 


FIG. 5. Panel (a): the diagrams for coupling between PDW and SC orders. Panel (b): the diagrams for coupling between 
CDW and SC orders. We have omitted many diagrams of the same type. Dashed lines stand for PDW order, wavy lines stand 
for CDW order, and double lines stand for SC order. 


Evaluating the integral in the same way as with did for Ii — I 4 earlier, we obtain [ 

(fk I 1 1 


h=-T 


=T 

= -T 


S/ 

E f dti 

J 


47r2 iwm - ei(fc) -iuim, - e 2 {k) [-iuim - e 5 (-fc )]2 
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(33) 


The second diagram in of Fig. [51(a) yields the coupling term that depends on the phases of -0 and T. This term is 
Combining contributions from the two types of diagrams, we obtain 

>Spdw/sc = 2|/4 |(i^cP + + 2\hUci>D{^*f+rcrD^% ( 34 ) 

Keeping in mind that \'ipc\ = \'^d\^ we find that 5pdw/sc is minimized when At the 

minimum, <Spdw/sc = 0? i-^-? PDW component of the mixed CDW/PDW order does not couple to SC order. 

For CDW components, the coupling terms are similar, as we show in Fig. [51(b). Following the same steps we obtain 

5cdw/sc = 2|/4 |(|m|" + + 2\h\{pAP*B+P*APBm^. (35) 

In distinction to PDW case, the relative phase of pA and is locked at ±7r/2 due to coupling to fermions away 
from hot spots Q. For such a phase difference, the second term in 5cdw/sc vanishes and the coupling between CDW 
component and SC reduces to 

Scd^/sc=Mh\{\pA\^ + \pBfm^. (36) 

Using Ip^I = \pb\ = |A|, we obtain 5i„t = >Spdw/sc + ^cdw/sc = /3'|A|2|^'|2 and /3' ~ \h\ ~ 1/ivjpT) ~ l/{vl.g). This 
is the main result of this analysis. 


II DETAILS OF THE FERMIONIC SPECTRA 

In this section we present the details of our analysis of the feedback from the CDW/PDW order on the fermionic 
dispersion in the antinodal region. 


A. CDW and PDW orders away from hot spots 

To calculate the fermionic dispersion in the whole antinodal region, we need to extend our analysis of the CDW and 
PDW orders to regions away from hot spots. We recall that in the SU(2)-symmetric hot spot model, CDW and PDW 
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FIG. 6. Schematic plot of the CDW and PDW magnitudes as a function of kx of momentum in the antinodal region around 
momentum (7r,0), obtained from Eqs. (l3^ and (H0|) at low temperatures T ei, 2 . Note that kx = tt — Ql2 is at the hot spot 
and kx — IS at the BZ boundary. 


components of the CDW/PDW order have equal magnitudes. Once we move away from hot spots, the equivalence 
gets lost. The momentum dependence of CDW and PDW order parameters p{k) = ^q,^(cJ^(Q/2 + k)cp{—Q/2 + /c)), 
and 'ip{k) = ia^^{ca{Q/2 + /c)c^(Q/2 — k)) at arbitrary k can, in principle, be obtained by solving integral equations 
in momentum and frequency. However, solving such an integral equation is a rather involved procedure. Here we 
adopt a simplified approach and just compare the kernels of the ladder equations for p{k) and 'ip{k) on both sides of 
a hot spot, i.e., near the zone boundary and closer to zone diagonals. 

Because in ladder series for p{k) and 'ip{k) each spin-fermion interaction changes fermionic momenta by tt = (tt, tt), 
we consider the effective kernel for small momentum transfer, made out of two subsequent terms in the ladder series. 
For CDW order, such term is 


Hp = g‘ 


j 

J UJirn ik 


G{ujfri’) Q/2 + TT + k + k)G{ujjji, —Q/2 + tt + k + f G{uj'^^ Q/2 + k + k'^Giuj'^^ —Q/2 + k + 

(37) 


where k is a deviation from a given k. We assume that relevant k are small enough and linearize fermionic dispersions 
in k. For PDW order, the product of two subsequent terms in the ladder series is 


— 9^ [ GiidjYi^ Q/2 + TT + k + k)G{uJrm Q/2 + TT + k — fc) [ Q/2 + k + k'^Gitu'^^ Q/2 + k — k')^ 

(38) 


The effective interaction g does indeed depend on momentum and also on frequency. However, when we consider 
fermions away from hot spots, this interaction is not singular and can be approximated by a constant. 

We first consider k at the Brillouin zone boundary, i.e., at k = (tt, 0) if Q = (0,(5). We linearize the fermionic 
dispersion as e(Q/2 + 7r + k + fc) = 62 + V 2 ky^ e(—Q/2 + 7r + k + ^) = €2 — V 2 ky^ e(Q/2 + k + fc) = ei — viky^ and 
e(—Q/2 -b k + ^) = ei viky. For antinodal fermions, we assume that 61^2 A. Note that the Fermi velocities 
between both fermion pairs that differ in momentum by Q are antiparallel. Evaluating the integrals, we obtain at 
small T : 
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FIG. 7. Panel (a): The dispersions of energy eigenstates in the presence of CDW and PDW orders for Domain II at kx = O.Qtt 
and p = = A = 80 meV. In general p and -0 are not equal and depend on momentum, but this does not change our conclusion. 

Panel (b): The spectral function in the presence of CDW and PDW orders for Domain II for the same parameters. Near the 
energy of the original band, the shadow bands and the original band strongly hybridize and the spectral weight of the original 
band is hence reduced. Panel (c): For comparison, the spectral function for Domain I at = O.Ott. Panel (d): The spectral 
function for Domain I and Domain II combined. The features from Domain I are clearly visible. 


Here A and A/, are energy and momentum cutoffs, respectively. We see that lip ^ H^. As the consequence, 
at the zone boundary the magnitude of the CDW order is much larger than for PDW order. This result can be 
straightforwardly understood because in the CDW channel Fermi velocities of the two fermions which form the 
condensate are antiparallel, while in the PDW channel they are parallel. 

As we move towards hot spots, the difference between the magnitudes of CDW and PDW orders gets weaker, 
and when k is between hot spots, CDW and PDW orders become degenerate. If we move further away from zone 
boundary, the kernels in PDW and CDW channels remain comparable to each other. Therefore, moving from the BZ 
boundary toward BZ center, the PDW gap should increase and finally become the same with CDW gap at the hot 
spot, as we sketch in Fig. [6l This is the main conclusion of this Subsection. 


B. Feedback from the mixed CDW/PDW order on the dispersion of antinodal fermions 


The feedback from a pure CDW order on the fermionic dispersion has been analyzed in Ref. 0. We follow the same 
approach detailed there, only in this work we also include PDW off-diagonal matrix elements. The matrix from whose 
determinant one can extract the fermionic energy contains the original band and the shadow bands with momentum 
shifted by ±(5,=1=2(5,.... For our purposes it is sufficient to only retain four shadow bands: two from CDW coupling 
with momentum (0, TQ), and two from PDW coupling with same momentum, because in the antinodal region all 
other shadow bands have high energy. Therefore, the energy matrix to diagonalize becomes 5x5 and we have 


E{k) = 


e{k) 

Phik) 
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0 
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(41) 


where Q = (0,(5). Since we are mainly interested in the fermionic spectra close to the Fermi surface, which is 
almost “horizontal” for the antinodal region close to ( 7 r, 0 ), we simply take p and ^ to be functions of kx^ i.e., 
PQik) = P-Qik) = p{k^) and -ipQik) = -tp-gik) = 'ip(k^). 

For the input values for CDW and PDW order parameters p{kx) and we use the results from the previous 

subsection. We set CDW order parameter p{k) to be independent on k in the antinodal region, and assume that 
PDW order parameter is zero at the Brillouin zone boundary and increases and approaches the CDW order as one 
moves the scan at fixed kx from the zone boundary towards the one which passes through a hot spot. Specifically, we 
set 


pikx) = A and ip{kx) = 


(42) 
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FIG. 8. The computed fermionic spectral function for kx = tt, kx = O.Ttt and kx = O.Gtt, assuming that only CDW order is 
present. The white parabolic line corresponds to the bare dispersion. Note that the dispersion crosses the Fermi level by a 
state moving from the upper branch, leaving a gap below the Fermi level. 


We take A = 80 meV to match the ARPES data. As experimental input, we use the dispersion in nearly opti¬ 
mally doped Pb-Bi2201 from Ref. 0: e{kx^ky) = —2t{coskx + cosky) — At {cos k^ cos ky) — 2t" {cos2kx + cos2ky) — 
At" (cos 2kx cos ky + cos kx cos 2 ky) — eo, with t = 0.22eV, t = —0.034315eV, t' = 0.035977eV, t" = —0.0071637eV , 
and eo = —0.240577eV. For this material, the CDW wave-vector is Q = 0.37r (see Ref. [ 1 ). We diagonalized the 5 by 
5 energy matrix (|Hj) and computed the spectral function /(cj, k) oc Im(( 0 /^( 0 ;) 0 ^( 0 ;))). Our results are shown in Fig. 
2 of the main text. 

Because the CDW/PDW order breaks C 4 symmetry, the measured spectra should be a combination of contributions 
from two domains Q. So far we have only considered one domain (Domain I) with the ordering momentum Q = (0, Q). 
For the other domain with ordering momentum Q = (Q,0) (Domain II) we found that for most of momentum range 
the CDW and PDW shadow bands and the original band have comparable energies [see Fig. [71[a)] in which case the 
spectrum is “dimmer” and features are much less pronounced. For comparison, in Fig. 0 we present the spectral 
function from kx = O.Qtt separately from the two domains and their combination. We see that the visible features 
are those from Domain I, i.e., the contribution from Domain II can be safely neglected. We found similar results for 
other values of kx. 


C. Incompatibility of the a pure CDW state with the ARPES data 

For comparison, we present the calculated fermionic spectra along multiple cuts like we did in Fig. 2(b) for a pure 
CDW state, and show that it is generally inconsistent with APPES data. Such an argument was first made by P. A. 
Lee 0, and we reproduce it here. 

We perform the same calculation as in previous Subsection, the only difference being that the PDW order parameter 
7 /; has been set to zero from the very beginning. We show the resulting fermionic dispersions for kx = tt, kx = 0.77r 
and kx = O.Gtt in Eig. [HI 

Comparing this with Eig. 2(c) in the main text, we see that at = tt our spectral function is generally consistent 
with experimental data. However, for other kx values, there is clear discrepancy. Namely, in the ARPES data, the 
Eermi arc region begins by “the lower branch moving up”, while in our result with purely CDW order, the Eermi arc 
begins by “the upper branch moving down”, leaving a finite gap below the Eermi level. This discrepancy is remedied 
by the presence of a PDW order, as we demonstrated in the previous subsection. 


[1] Y. Wang, D. F. Agterberg, and A. V. Chubukov, Phys. Rev. B 91, 115103 (2015). 

[2] D. Chowdhury and S. Sachdev, Phys. Rev. B 90, 134516 (2014). 

[3] Y. Wang and A. Chubukov, Phys. Rev. B 90, 035149 (2014). 

[4] K. B. Efetov, H. Meier and C. Pepin, Nat. Phys. 9, 442 (2013). 

[5] M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075128 (2010). 







17 


[6] L. E. Hayward, D. G. Hawthorn, R. G. Melko, and S. Sachdev, Science 343, 1336 (2014). 

[7] Rui-Hua He, et ah. Science 331, 1579 (2011). 

8 W.D. Wise, M.G. Boyer, K. Ghatterjee, T. Kondo, T. Takeuchi, H. Ikuta, Y. Wang, E.W. Hudson, Nat. Phys. 4, 696 (2008). 
[9] P. A. Lee, Phys. Rev. X, 4, 031017 (2014). 


